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ABSTRACT 

We analyse the full shock formation process in electron-ion plasmas in theory and simulations. 
It is accepted that electromagnetic shocks in initially unmagnetised relativistic plasmas are trig¬ 
gered by the filamentation instability. However, the transition from the first unstable phase to the 
quasi-steady shock is still missing. We derive a theoretical model for the shock formation time, 
taking into account the filament merging in the non-linear phase of the filamentation instability. 
This process is much slower than in electron-positron pair shocks, so that the shock formation is 
longer by a factor proportional to y/rriijrnehi{mi/me)■ 

Subject headings: instabilities — relativistic processes — shock waves 


1. Introduction 

Collisionless shocks are ubiquitous in astro- 
physical environments, such as gamma-ray bursts, 
active galactic nuclei or pulsar wind nebulae. 
They are especially important in the context 
of cosmic ray acceleration. Once the shock is 
in a quasi-steady state, the jump conditions 
can be determined from the conservation of 
mass, energy and moment um in a fluid model 
( Blandford fc McKed 1976 1. The density jump 
from the upstream to the downstream of a rela¬ 
tivistic strong shock is given in 2D by n 2 lni ~ 3. 
Plasma instabilit ies are the med iators of such col¬ 
lisionless shocks (Sagdeey 1966ll . but it is still not 


^also at DCTI, ISCTE - Lisbon University Institute Por¬ 
tugal 


known precisely how the transition from the initial 
plasma turbulence to the final quasi-steady state 
of the shock happens. 

In a symmetric counter-streaming flow of 
charged particles, plasma instabilities develop, 
causing an isotropisation of the particle momenta 
and initiating the deceleration of the plasma flows. 
In relativistic, cold and initially unmagnetised 
plasmas, electromagnetic current filamentation 
modes are dominant w hich inhabit a strong per - 
pendicular component (Weibel I959I : Friedllig^ . 
The energy is transferred from an initial longitu¬ 
dinal streaming into the perpendicular directions. 
The particle forward motion is slowed-down and 
a collisionless shock starts to form with a density 
ratio of n^jn-i > 2. 

Particle-in-cell (PIC) simulations are an ex- 
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cellent tool to investigate the shock formation 
process, since non-linear processes are involved 


process, since non-linear processes are involved 

(iMedvedev et al ■I 2 OO 5 I: iHededal et aI1l2008l: ISpitkovskvl 


the analytical evaluation of the shock formation 
time. As already highlighte d by various authors 
(iLviibarskv fc Eichleil 20061 : Shaisultanov et al 


Fiu^ae^_alJl2012 ^ Sironi^_^ ■ | 2013t Stockem et al. 


200 J Martins_^talJ_ 2009^ Nishifewc^^^jb[ 20nl I 2 OI 2 I : Davis et al.ll2013ll , the electron Weibel insta- 


2014 iHuntington et al.l l2015ll . In a previous 
study, we identified the shock formation time of 
electron-positron pair shocks as twice the satu¬ 
ration time of the magnetic field amplification 
due to the filamentati on instability, T/,e = ‘^Tg^e 
( Bret et al. 20131 2014ll . Surprisingly, electron-ion 
shocks do not show the same feature as one would 
expect from the rapid relativistic mass increase 
of the electrons. Our analysis shows a shock for¬ 
mation time Tf^i Ri 3 However, a 

similarity between both scenarios was observed: 
the steady state shock is formed when the first 
ions start to recirculate. The non-linear phase is 
faster in pair shocks and no additional time is nec¬ 
essary for the merging of magnetic flux tubes, so 
that the recirculation process starts earlier. As a 
consequence, particles are mainly scattered in the 
turbulence rather than doing a full gyro rotation, 
which is why the isotropisation process is less ef¬ 
fective. We present a detailed theory for the shock 
formation in electron-ion plasmas, and compare it 
to state-of-the-art PIC simulations. 

2. Shock formation time 

We consider a simple scenario for shock for¬ 
mation with the plasma initially being unmagne¬ 
tised and symmetric counterstreaming relativistic 
beams of electrons and ions of mass ratio > me 
and Lorentz factor 70 . We refer to pair shocks 
of electrons and positrons when mi = me- The 
beams are initially cold, characterised by the tem¬ 
perature parameter /x = /ksT, where the 

mass m and temperature T refer to the respec¬ 
tive species. Such scenarios are a stimulating en- 
viro nment for t he current filamentation instabil¬ 
ity ( Fried 19591) to occur. The growth rate of 
the cold current filamentation instability, for elec¬ 
tron beams as well as for ion beams, is given by 
(5 = \/ 2 / 7 o/ 3 qWp, where /3o = vojc denotes the nor¬ 
malised fluid velocity of the beam and the plasma 
frequency is given by ujp = ^/ 47 ^nne^/m with th e 
initial uniform beam density uq { Bret et al. 2010l) . 

Let us start describing the mechanisms at work 
that lead to shock formation, before we turn to 


bility is first triggered when the two plasmas start 
overlapping. By the time it saturates, it has gen¬ 
erated filaments of the size of the electron Larmor 


radius in the field at saturation (jBret et al.ll201,‘ll) . 
At the time the ion Weibel instability starts to 
grow, an unstable wavelength has already been 
seeded. It corresponds precisely to the one that 
has resulted from the merging of the filaments as¬ 
sociated with the electron Weibel instability with 
a typical length scale of c/wpe. As a consequence, 
by the time the ion Weibel instability saturates, 
the field is near equipartition with the ions, but 
the filaments are still the size of the electronic Lar¬ 
mor radius, not the ions one, since the instability 
on the ions was initially seeded at these length 
scales. Such filaments are too small to efficiently 
deflect the ion flow, an d need to merge in order 
to reach required size jMilosavlievic et al. _ 200d 


2nnsl: 


Davis et al 


Spitkovskvl I2008I iChang et al 
2 OI 3 I ). Once the appropriate number of merging 


events has been achieved, the filaments reach the 
size of the ion Larmor radius. Only then is the 
ion flow deflected enough for the shock to start 
forming. Assuming this happens at a time r and 


follow ing the reasoning explained in I Bret et al 


( 2014 ). the shock formation time will be 2 t in 2D, 
and 3r in 3D. 

2.1. Saturation phase 


The electron Weibel instability grows first, am¬ 
plifying the field from its fluctuation value up to 
nearl y equipartition with the electrons (I Bret et al 
I 2 OI 3 I ). The saturation time is here given by 


= Hw, 


pe ; 


( 1 ) 


where H is the number of e-foldings of the insta¬ 
bility and SeUJpe growth rate. By this time, 
the field has grown to nearly equipartitio n with 
( Medvedev &: Loe^ll999t Isilva et al. 2003 ) 


g = STrynomeC^ 


( 2 ) 


From the cyclotron frequency = QBg^e/j'meC, 

we derive the size of the filaments at saturation, 
which is also the electronic Larmor radius in B, « 
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(we set vq ~ c) 




7 c 
2 


(3) 


A full quantitative understanding of the saturation 
magnetic f ield is still lacking and will be explored 
elsewhere ( Schoeffler et alJ[2015ll . 

Contrary to the electron current filamentation 
instability, which has to amplify the plasma ther¬ 
mal fluctuations, the ion Weibel instability finds 
an unstable mode already seeded, and further am- 
plihes it, growing preferentially at that wavenum¬ 
ber. This is clear from the fact that the growth 
rate of the ion Weibel instability depends slowly 

w r /I 

on k for wavenumbers k > 7- 1 -|- 

( Davidsonlll972 ). At saturation, the ion Weibel 
instability has grown the field to 


l|2005l) who considered a simple 2D model of in- 
hnitely long cylindrical filaments of radius D/2, 
spaced by the distance i^l. 

The filaments need to merge n times in order to 
grow from the electronic Larmor radius (l3l) to th e 


ionic one. According to iMedvedev et al.l (1200511 , 


the expression for the merging time is different 
whether the transverse motion is relativistic or not 
(merging of the filaments implies a motion trans¬ 
verse to the flow). At any rate, the maximum 
transv erse velocity achieved during one merging 


reads ( Medvedev et al. 2005l l 


'^m 


1.67 Dujpi 
4 ’ 


(7) 


considering vq ^ c. We now compare this ve¬ 
locity to c. For the merging motion to be non- 
relativistic, we would need. 


= STrynom^c^, (4) 

while the size of the filaments is still given by Eq. 
©. The growth time of the ion phase of the in¬ 
stability is 



(5) 

The total duration of this saturation phase is Ts = 
Ts.e+Ts.i, that with S = 6e = Siy/rrii/me we obtain 

+ ljw7- (6) 

Considering a number of e-foldings 11 ^ 20, we 
hnd that even for an artihcially low mass ratio of 
100 the term 11mg/mi/ In ^Jrai/rrie 1. As 
a consequence, the duration of the whole satura¬ 
tion phase simply reads ^ Tg^i. In reality, the 
ion growth rate 5i is lower than the 

growth rate of the cold ion current filamentation 
instability, because it grows over a background of 
hot electrons. Yet, the consequences are negligi¬ 
ble, see section [231 



2.2. Filaments merging phase 

The dynamics of the merging 
has been studied previously in 


of the filaments 
Medvedev et al. 


Vm c ^ D <s: 2 . 4 Y 7 —. ( 8 ) 

^pi 


If this condition is fulhlled at the end of the merg¬ 
ing phase, then it is always fulfilled, since filaments 
are growing with time. Replacing therefore D by 
the ion Larmor radius Li = y^'y/2c/u}pi, we find 
the condition above is always satisfied. 


For t he regime of no n -relat ivistic transverse 
merging, Medvedev et al.l ( 2005ll established that 
all merging events are taking approximately the 
same time 

ro = 23/2^a;-4. (9) 


This number now has to be multiplied by the 
total number n of merging events. Since the 
initial and final filaments size are the electrons 
and ions Larmor radii respectively, n is given by 


( Medvedev et al. 2005ll 


Li = 2”/2 r ^ n = 2 


\n{Li/Le) 
In 2 


ln(m,/me) 
In 2 


( 10 ) 

The total duration of the merging phase is there¬ 
fore given by 


23/2 

Tm = riTo = ln(mi/me)7^/^w“4. (11) 


^The notation in this p aper differs from the one used in 
IMedvedev et all ll 2005 h . The correspond ence is obtained 
replac ing the quantity Iq / y/JiQ, found in IMedvedev et all 
j 2005 t) , by Dijpivo/Ay^. 
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2.3. Shock formation time 

We can now proceed to the calculation of the 
full shock formation time. As soon as the filaments 
are large enough, the ions are efficiently deflected 
by the field over a distance close to the ion Larmor 
radius. This happe ns at a time r,„ + Tm. ■ Assuming 
then = \/jj2 (jBret et al. 20ir)ll , we find, 


H" '^m — 


= {. 




V23/2 ' in2 J 

4.437'/^ In 


-OJ. 


pi 


( 12 ) 


pair p lasma by setting me = mi ( Bret et al ] |2013L 


pair p i; 

l2ni4l l. The reason for this is that here we simply 
neglect the electron instability phase in the end re¬ 
sult. Finally, we note that shocks in electron-ion 
plasmas form much slower than in pair plasmas 
for two reasons: on the one hand, the instabil¬ 
ity mechanism is slower while on the other hand, 
the merging phase is negligible in pair plasmas be¬ 
cause the instability already generates filaments 
and fields large enough to deflect the flow. 

2.4. Comparison with the pair shock for¬ 
mation time 


We thus find that the merging time is longer 
than the ion Weibel saturation time by a factor 
2^/ In 2 ^ 11.5. As a consequence, the uncertainty 
on the ion Weibel growth rate does not affect sig¬ 
nificantly the formation time. In order to find out 
whether or not the flow is stopped in the over¬ 
lapping region by -I- we need to compare 
L = 2c{Ts+Tm), the size of this region at this time, 
with the ion Larmor radius Li. It is straightfor¬ 
ward to show that 

L/Li = 8.83^2 In—» 1. (13) 

me 


As a consequence, the incoming flow is stopped 
in the overlapping region, and the shock starts 
forming. The downstream density at time Tg + Tm 
is still only twice the upstream density. In or¬ 
der to reach the expected density jump of 3 for 
the 2D cas^l, the system needs to evolve another 
Ts + Tm since the overlap ping region no longer 
expands (|Bret et al. 2014 ). In the 3D case, the 
expected density jump is ~ 4 so that waiting 
another 2 (ts + Tm) is required to bring enough 
material in the central region. Note that the 
present 2D model of the filaments merging as 
been successfully tested ag ainst 3D PIC simula¬ 
tion ( Medvedev et al.ll2005r) . The shock formation 
time Tf^i hnally reads. 


iLi 

d 


4.437^/^ In 


me 


(14) 


where d is the dimensionality of the system (d = 
2(3) for a 2D (3D) setup). Noteworthily, the ex¬ 
pression does not reduce to the one obtained for a 


^ See IStockem et al.l 120121 ) and I Bred 12014 ) for the validity 
of the Rankine-Hugoniot jump conditions in collisionless 
plasmas. 


For comparison of the shock formation times 
of pair and electron-ion shocks, Eq. (fTT) gives the 
form ation time Tj^e for the pair case (jBret et al. 
20141) 

Tf^e = dTs,e = d5e^IlLO~e (15) 

in terms of the dimension d of the problem, 
whereas from Eq. 0 we obtain for electron-ion 
shocks 

^ ln{mijme)\/mijme. (16) 

T-f,e n 

Considering 11 ^ 12, which is the value obtained 
in the pair shock simulations presented in the next 
section, we find that an electron-ion shock forms 
60 times slower than a pair shock for a mass ra¬ 
tio of 400, and 166 slower for a realistic mass ra¬ 
tio. For a mass ratio of 400, Tf^i = Si/miJm^Tf^e 
whereas = 2.4y^mi/meT/,e for a mass ratio of 
100 . 

3. Discussion of the results 

We present now a series of 2D PIC simula¬ 
tions in order to validate the theoretical model. 
The counterstreaming beams were simulated in 
a simulation box with perfectly reflecting wall 
in the longitudinal direction and periodic bound¬ 
ary conditions transversally. The bulk is prop¬ 
agating along the xi axis with Lorentz factors 
7o = 25—103, mass ratios mi/me = 50 — 400 
and = 10®7o. The simulation box dimensions 
are = 450^/^c/a;pe and Ly = 150 y^c/tUpe 
with Ax = Ay = 0.05 y/^c/ojpe- 

The magnetic field energy density is plotted in 
Fig. [T] in order to get information about the role 
of the filamentation instability as mediator of the 
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shock formation. This was done for a small slab 
along xi in a region close to the wall, which will 
later be a region in the far downstream. 



Fig. 1.— Normalised magnetic energy density for 
an electron-ion shock with Wi/me = 400 (blue) 
and a pair shock (black) with 70 = 25. The lines 
indicate the saturation times of the filamentation 
instability and the shock formation times. 


In the case of the pair shock the linear phase 
of the instability - where the magnetic field grows 
exponentially in time and the magnetic field en¬ 
ergy density is oc exp(2^ei) “ can be clearly 
distinguished from the non-linear phase where the 
magnetic field has saturated. The growth rate of 
the cold electron instability fits very well the the¬ 
oretical value Se = 0.28 ujpe, as well as the satura¬ 
tion field Bf^e ~ 7 meCUJpe/e. The predicted satu¬ 
ration time of the filamentation instability of the 
pair shock Ts,e = 40a;“g^ and steady state shock 
formation time ry_e = 80 uj~^ are shown in Fig. [T] 
and match well with the simulation data. 


In the case of the electron-ion shock with 
rriilme = 400, the evolution of the magnetic en¬ 
ergy density shows several stages. The ion Weibel 
instability grows slower than 5i = yra^Jrriide = 
0.014 Wpe du e to the influence o f the hot electron 
background (jShukla et al.l 1201211 , which is negli¬ 
gible when compared to the full formation time. 
The final magnetic field at saturation of the ion 
instability matches approximately the theoretical 
value Bf^i Ri 140cme<jjpe/e as well as the satura¬ 
tion time in the simulation Ts^i = 18a;“4. The 
theoretical model predicts a saturation time of 
13 Wp/, see Eq. ([51). 


The black lines in Fig. [5] indicate the shock 
front of the steady state shock with a jump of 
= 3. This line is extrapolated to cci = 0 
in order to define the shock formation time rj as 
the intersection with the time axis. For the pair 


shock (Fig. [2^) it matches the theoretical value 
= 2Ts_e) while for the electron-ion shock we 
observe = 226 > 2x^7 = 36u;“4 (Fig. 

ISb)- A systematic study with different velocity 
and mass ratio parameters provided a factor of 
2.5 mil me with respect to the shock formation 
of pair shocks, which is consistent with Eq. (HI 
(see Fig. El). The shock width imposes an uncer¬ 
tainty in the determination of the shock formation 
time, which is represented by the error bars in Fig. 

El 




lo 


Fig. 2.— Ion density against space Xi and time 
t in units of w”/ = ^rmlme for 70 = 25, 
rriilme = 1 (a) and mijme = 400 (b). The black 
line indicates the shock front with n^lrix = 3 and 
the pair shock for comparison (dashed blue). 
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Fig. 3.— Shock formation time v s. 7n for a 


2D p air shock with mi/me = 1 (black) (|Bret et al 


20141 ) and electron-ion shock with different mass 


ratios (red). The error bars determine the uncer¬ 
tainty due to the finite size of the shock front. 


The delayed shock formation process in electron- 
ion shocks due to the merging of the filaments is 
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demonstrated in Figs. |4] and [5] In both cases, the 
accumulation of particles becomes very effective 
at the time when the phase space of the different 
beams starts to mix (Fig. U). The first ions start 
to recirculate and change the sign of their longi¬ 
tudinal momentum. In panels (c) and (d) of Fig. 
m this is demonstrated by a small fraction of par¬ 
ticles close to zero momentum. At this stage the 
magnetic field turbulent scales are large enough 
that the ions can finish at least half a gyro circle, 
and the accumulation of particles becomes effi¬ 
cient. In pair shocks this process happens already 
right after the saturation of the filamentation in¬ 
stability, at tuipe = 65, whereas in electron-ion 
shocks this process takes much longer. For a mass 
ratio of milrrie = 400, the efficient gyro reflection 
was observed only at tiOpi = 226, which is slightly 
before the steady-state shock has formed. 




[c/mj 
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Fig. 4.— Momentum phase spaces of a pair 
shock at 65 Wp/ (a) and electron-ion shock with 
nii/me = 400 at 226 w”/ (b) with respective aver¬ 
aged particle distributions (c) and (d). 


An analysis of the magnetic field structure con¬ 
firms this result. At the time when the magnetic 
field saturates, the transverse size of the mag¬ 
netic field flux tubes is still too small in order 
to scatter the ions efficiently. The particles will 
feel the impact of another flux tube long before 
they can finish a gyro circle. For pair plasmas 
(Fig.[S^) the magnetic filaments at saturation time 
Ts,e = 40a;have reached a transverse spread of 4 
c/a;pe, while the maximum magnetic field strength 
is of the order of = ^meCLOpeje. Particles with 
Lorentz factors 7 < 12 have Larmor radii that fit 
into this scale, meaning that they can recirculate 
in the field before being deflected by a different 
flux tube. This situation is different for electron- 


ion plasmas. For comparison, we plotted the mag¬ 
netic field structure at Tg^i = 18 w”/ in Fig. [Dd. 
At this stage, the magnetic field filaments show 
a transverse size on the order of 2 cjujpi while 
B 3 Ri 5meCU}pe/e, so that only particles with non- 
relativistic Lorentz factors 7 ~ I will have Larmor 
radii on the same scale as the magnetic field tubes. 
The filaments have to undergo a further merging 
process as described in Eq. (HU until the trans¬ 
verse filament size becomes of the size of the ion 
Larmor radius for particles with 7 > 25. 



X|[i:/l0„l X, [c;m„| 

Fig. 5.— Magnetic field of the pair shock at satu¬ 
ration time Ts^e = 40 0 ;“^^ (a) and for electron-ion 
shock at Ts,i = 18 a;”/ (b). 


4. Conclusions 


We have investigated the full shock formation 
process in electron-ion plasmas in theory and sim¬ 
ulations. In contrast to electron-positron pair 
shocks, where the shock formation time was found 
to be twice the saturation time of the filamen¬ 
tation instability, = 2 Ts^e^ for electron-ion 
shocks the process is delayed to approximately 
Tf^i = The shock formation time 

is a sum of the saturation time of the instability 
plus an additional merging time, coming from the 
merging of filaments to the ion Larmor radius. At 
the time of the saturation of the instability, the 
filament size is still on the order of the electron 
Larmor radius. An extra time is thus necessary 
for the condition for shock formation to be met, 
which is not t he case for pair shocks. We applied 
the theory bv iMedvedev et ah ( 2005 1 to predict 
the merging time r™. The merging time retrieved 
from 2D PIC simulations is in agreement with the¬ 
ory and the shock formation time was confirmed to 
be 2(Tg^i + Tjn). Slightly before this time, the first 
recirculation of ions was observed. At this stage, 
the scale of magnetic turbulence is large enough to 
trigger the density compression that precedes the 
full shock formation. 

The different scales of the magnetic turbulence 
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in electron-ion shocks compared with pair shocks, 
might have consequences for the particle accelera¬ 
tion process. This will be investigated in our fol¬ 
lowing project. 
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